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TEV BLAZARS AS SEEN BY THE CAT TELESCOPE 

F. Piron 1 and the CAT collaboration 

Abstract 



To date, only two extragalactic objects have been firmly established as 
very high-energy gamma-ray sources in the Northern sky: these are the two 
blazars Markarian 501 and Markarian 421. This paper reviews the most stri- 
king results obtained from these sources by the CAT atmospheric Cherenkov 
imaging telescope, with a particular emphasis on the 1999-2000 and 2000- 
2001 observation campaigns of Markarian 421. 



<^> 1 Introduction 
Oh 



Blazars are believed to be radio-loud active galactic nuclei (AGN) which have a 
relativistic jet with its axis closely aligned to the line of sight. Due to Dopplcr 
boosting, the jet emission dominates that of the accretion disk (and of the host 
galaxy) over a large energy domain. It can extend up to the very high-energy 
(VHE) range (>100GeV) in some cases. TeV blazar observations offer the pos- 
sibility of investigating the physics of jets more deeply, including particle accele- 
ration and energy extraction in the vicinity of the AGN's central "engine" . The 
identification of the VHE 7-ray production region with the base of the jet is indeed 
suggested by their highly variable emission, which has been already observed down 
to sub-hour time-scales [0, thus indicating a small spatial extension via a simple 
causality argument. 

Since the detection by the Whipple Observatory of Markarian 421 (Mkn 421) 
]l"l| and Markarian 501 (Mkn 501) Q as the first extragalactic TeV sources, 
blazars have been intensively observed by ground-based Cherenkov atmospheric 
7-ray detectors. The CAT (Cherenkov Array at Themis) imaging telescope Q, 
which operates above 250 GeV since Autumn 1996, has studied Mkn 421 and 
Mkn 501 in different states of activity. The complete analyses are detailed in the 
literature || [l0|, ^| ; here we present the main results of this work. 
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Fig. 1. (a) Mkn 501 and Mkn 421 nightly-averaged integral flux above 250 GeV as 
recorded by CAT during their periods of most intense activity. Dashed lines show the 
mean flux for each observation period; (b) Zoom on the nights between 3 and 6 February, 
2000. Dashed lines show the Crab nebula flux; (c) VHE spectral energy distributions 
(SEDs) of Mkn 501 (hatched areas, multiplied by 10 for clarity) and Mkn 421 (unfilled 
areas). Each area stands for the 68% confidence level contour (see [JlOj for details). 



2 VHE properties of Mkn 501 and Mkn 421 

The light curves of these two blazars, as sampled by CAT during their most active 
periods, are shown in Fig. 0(a). Both sources are highly variable and a few series 
of bursts can be distinguishedR For instance, the flux of Mkn 421 was measured 
at a level of (67.3 ±2.6) x 10""cm" 2 s _1 above 250 GeV during the night of March 
23/24, 2001 (JD 2451992): this is comparable to the highest TeV flux ever recorded 
by CAT, i.e., that of Mkn 501 in the night of April 15/16, 1997 (JD 2450554), 
which reached (83.2 ± 3.8) x 10~ n cm -2 s _1 . During these flaring periods both 
sources underwent a large night-to-night variability on several occasions. Mkn 421 
is, however, unique since it has also exhibited important flux variations within a 
single night. As an example, its light curves for three nights from the 3 rd to the 5 th 
February 2000 are shown in Fig. |l|(b). While the fluxes recorded by CAT during 
the first and last nights were stable, the source activity changed dramatically in a 
few hours during the second night: it was reduced by a factor of 2 in 1 hour and 
by a factor of 5.5 in 3 hours. Consequently, the 7-ray emitting zone must be very 
compact, with a size <10 light-hours if one assumes a typical value of 10 for the 
geometric Doppler factor (which reduces the time-scale in the observer frame). 



The uneven sampling of Mkn 421 light curve in 2001 is due to bad weather conditions. 



F. Piron : TeV blazars as seen by the CAT telescope 
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Fig. 2. (a) Spectral measurements on Mkn 421 by HEGRA, CAT and Whipple between 
1995 and 2001. We added the CAT results obtained during the nights of February 4/5, 
2000, and March 24/25, 2001. The results are given in the {<f>o,y} plane of spectral 
parameters in the power-law hypothesis. All error bars take account of both statistical 
and systematic errors (see jl^] for details and references therein). The xV n corresponding 
to the absence of any dependence of the spectral index with flux is 3.5; (b) Mkn 501 X-ray 
and VHE SEDs for April 7 th and 16 th , 1997, as simultaneously measured by Beppo-SAX 
and CAT Ji], . Full lines come from our homogeneous SSC model |J . 

As shown in Fig. |l|(c), Mkn 501 clearly showed a curved spectrum in 1997, 
with a 7-ray peak lying above the CAT threshold. Mkn 421 is less extreme: its 
1998 time-averaged spectrum can be well fitted by a simple power-law with a 
differential index 7=— 2.88 ± 0.12 stat ± 0.06 syst , but there is some indication of 
curvature for the 2000 spectrum due to a larger signal-to-background ratio [|nj. 
The even higher statistics gathered in the sing le night of March 23 rd , 2001, yielded 

a more accurate determination of this curvature, with a differential spectrum — — oc 

dE 

U1 -2.77±0.09 stat ±0.06 syst -(0.72±0.27 3tat ±0.03 syat ) log ln BtcV ca-11 r • tt i i- 

E TcY K ' S1 ° . Still preliminary at the time 

of the conference, this result has been confirmed since then §. 

The question of spectral variability of TeV blazars has been the subject of a long 
debate. The hardening of the Mkn 501 spectrum during flares was firstly discovered 
by CAT on the basis of a hardness ratio during its 1997 flaring period It is 
illustrated in Fig. g(c) by the shift in the peak energy of the SED: 0.56±0.10TeV 
for the time-averaged spectrum and 0.94±0.16TeV for the highest flare of April 
16 th . Recently, the HEGRA experiment observed a similar feature by comparing 
the 1997 and 1998-99 time-averaged spectra jl). To address this question concer- 
ning Mkn 421, we compiled its spectral measurements in Fig. ||(a), in the {4>o,j} 

plane of spectral parameters assuming a power-law shape — — = ^o-^TeV (this 
hypothesis is justified even in the case of a curved sectrum, see ]1C|] for details). 
Despite the general trend seen in this figure, definite conclusions on a possible 
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spectral variability would be premature: simultaneous observations of more flares 
by different Cherenkov telescopes are still needed in order to exclude possible 
systematic effects between experiments (see Q for a discussion). 

3 Conclusion 

The multi-wavelength SED of Mkn 501 can be explained by leptonic models 0. 
The best example to date is shown in Fig. ||(b) , where a homogeneous Synchrotron 
Self-Compton (SSC) model fits well the X-ray and VHE spectra obtained by the 
Beppo-SAX satellite and CAT for two dates in April 1997 ||. In this framework, 
the CAT observations of Mkn 421 indicate that the peak energy of the inverse 
Compton contribution to its SED is significantly lower than the telescope detection 
threshold. This is not surprising since (i) the corresponding synchrotron peak is 
known to be lower than that of Mkn 501, and (ii) leptonic models predict a strong 
correlation between X-rays and 7-rays. This feature is illustrated in Fig. ||(b) by 
the correlated spectral hardening in both energy ranges between the two dates. 
As for Mkn 421, a correlation has been observed many times 0, [IJJ, but always 
in terms of integrated (and not differential) fluxes due to the lack of statistics. In 
addition the VHE spectral variability of Mkn 421 has not been clearly proven yet. 

Moreover, alternative scenarios based on an e _ p plasma and proton-induced 
cascades in the jet are still successful in interpreting TeV blazars' SEDs [jl3| |) . The 
study of the dynamic aspects of blazar jet emission, including the temporal and 
spectral correlations between various wavelengths, is thus required in the future 
in order to accurately constrain existing models and to understand the particle 
acceleration and cooling processes occurring at the sub-parsec scale in jets. It 
should help in discriminating between these models and allow, in particular, to 
address more deeply the crucial problem of the plasma jet content. 

References 

Aharonian, F.A., et al. 2001, ApJ 546, 898 

Barrau, A., et al. 1998, NIM A416, 278 

Djannati-Atai, A., et al. 1999, A&A 350, 17 

Gaidos, J.A., et al. 1996, Nature 383, 319 

Ghisellini, G., et al. 1998, MNRAS 301, 451 

Khelifi, B., et al. 2001, 27 th ICRC (Hamburg) 7, 6001 

Maraschi, L., et al. 1999, ApJ 526, L81 

Muecke, A., & Protheroe, R.J. 2001, APh 15, 121 

Piron, F. 2000, Ph.D. Thesis, University of Orsay-Paris XI 

Piron, F., et al. 2001, to be published in A&A 

Punch, M., et al. 1992, Nature 358, 477 

Quinn, J., et al. 1996, ApJL 456, 83 

Rachen, J. P. 1999, GeV-TeV Gamma-Ray Astrophysics Workshop (Snowbird), 41 
Takahashi, T., et al. 1999, APh 11, 177 



